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We explicitly construct a monopole creation operator: its vacuum expectation 
value is an order parameter for dual superconductivity, in that, if different from zero, 
it signals a spontaneous breaking of the U(l) symmetry corresponding to monopole 
charge conservation. 

This operator is tested by numerical simulations in compact U(l) gauge theory. 
Our construction provides a general recipe for detection of the condensation of any 
topological soliton. In particular our operator can be used to detect dual supercon- 
ductivity of the QCD vacuum. 
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I. INTRODUCTION 



A possible mechanism of colour confinement in Quantum Chromo-Dynamics (QCD) is 
dual superconductivity of the vacuum 

According to this scenario, the chromoelectric field is channelled into Abrisokov flux 
tubes ||, in the same way as the ordinary magnetic field is in superconductors: the word 
dual indicates the interchange of roles between electric and magnetic fields and charges. The 
chromoelectric field mediating the force between coloured particles is squeezed by Meissner 
effect into flux tubes of constant energy per unit length, giving rise to the confining linear 
potential. These flux tubes behave as strings The existence of strings in hadronic 

physics is supported by phenomenology |7]f|]. They have also been visualized by numerical 
simulations of QCD on the lattice JD|JT0[] . Some evidence in favour of dual superconductivity 



of the vacuum has been produced by Montecarlo simulations [11 



However a clear cut test of the mechanism is still missing. Such a test would be the 
detection of monopole condensation in the ground state, analogous to the condensation of 
Cooper pairs in the ground state of an ordinary superconductor. Condensation implies that 
the vacuum is a superposition of states with different charge, which, in turn, is nothing but 
a spontaneous breaking of the U(l) symmetry related to charge conservation Such 



a breaking is signaled by a non- vanishing vacuum expectation value (vev ) of any opera- 
tor carrying non-trivial charge, as, e.g. the scalar field of the Landau- Ginzburg model of 
superconductivity. That vev is called a disorder parameter in the language of statistical 
mechanics. 

In this paper, we give an explicit form of the creation operator of monopoles, together 
with a lattice version of it, in view of using its vev as a disorder parameter monitoring dual 
superconductivity. The operator is checked in U(l) compact pure gauge theory on a lattice, 
where the behaviour of monopoles is known p^JT^] . 

For QCD, the strategy is the following. The monopoles which are expected to con- 
dense and generate superconductivity are Dirac monopoles of a residual U(l) symmetry, 
which survives after a suitable gauge fixing, known as abelian projection 0. An abelian 
projection is defined as the gauge transformation which diagonalizes any operator trans- 
forming in the adjoint representation of the gauge group. The points where two eigenvalues 
of that operator coincide correspond to singularities of the gauge transformation and have 
the topology of world-lines of point monopoles fl,15l . Such monopoles have been observed 



on the lattice []T6|-20l. Of course the location and the number of monopoles do depend on 
the choice of the operator used to define the abelian projection: a possibility, advocated by 
't Hooft, is that physics, i.e. monopole condensation, is independent of it. The relevant 
abelian degrees of freedom can also be fixed by a somewhat different procedure, known as 



maximal abelian projection [fLq| . Our operator will allow to investigate unambigously what 
abelian projection, if any, defines the monopoles relevant to confinement. 

In this paper, we will only present the construction of the operator and the way it works 
in lattice U(l) compact gauge theory. For that theory, there exists a construction of a 
disorder variable describing monopole condensation |[21|| , which is rigorous but based on a 
particular form of the action (the Villain action). Our operator coincides with the above one 
in the Villain case, but can be used with different forms of the action. This is particularly 
important in order to use it for non-abelian theories, where the effective U(l) action after 
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abelian projection is not known. The study of monopole condensation in QCD is in progress 
and will be presented elsewhere. 



II. THE MONOPOLE CREATION OPERATOR 

For the sake of definiteness, we shall consider here only the case of a U(l) monopole. 
Our procedure can easily be generalized to all kind of solitons. 

Let £? M (x;y) = (0, b(x, y)) be the classical field produced in the location x by a Dirac 
monopole at rest in y. We can make any choice for the gauge, e.g. by putting the string 
along the positive z-axis. Then, defining r = x — y and r = ||r||, we have: 

T ' 

Mx, y) = g — - — r (i) 

r(r — r 3 ) 

where g is the charge of the monopole, satisfying the Dirac quantization condition: 

eg——, n integer (2) 

If Ilj(x, t) is the conjugate momentum to Aj(x, t), then the operator 

fj,(y, t) = exp{i J d 3 x ^(x, y) IL.(x, £)} (3) 

creates a monopole at the location y and time t. This can be immediately seen in the 
Schrodinger representation of the fields, where we have: 

My, 0) | A(x, 0)) = | A(x, 0) + b(x, y)) (4) 

Equation is a trivial consequence of the canonical commutation relations between the 
fields and their conjugate momenta, and is nothing but the field-theoretic equivalent of the 
familiar statement that e ipa translates the coordinate q by a: 

e ipa \q) = \q + a) (5) 

//(y, t), applied to any field configuration, adds a monopole to it. 
This can be restated in terms of commutation relations as: 

[Ai(x, t),n(y, t)] = k(x, y)n(y, t) (6) 
[n i (x,t),/i(y,t)] =0 (7) 

which also show that the electric field of the configuration is left unchanged. It is worthwhile 
to notice that the specific choice of the gauge for b is irrelevant: what really matters here 
is topology, which is independent of it. Our operator [i is similar to operators introduced in 
the literature by different constructions, in various contexts |22| . 

Now, if the ground state of the theory has a definite monopole number N, then, under a 
magnetic U(l) rotation, we have: 
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U\0) = e i(fN \0) 



Hence: 



which implies: 



Therefore if 



<0|//|0> = (0|/i|0)e 



Ho> = o 



up 



Ho>^o 



(8) 
(9) 



(10) 



(11) 



(12) 



then |0) is not U(l) invariant and there is spontaneous symmetry breaking of U(l). A 
translation by a static field g(x) such that curl g = in this language corresponds to a pure 
gauge transformation: 



7(t) = exp ji J d 3 y &(x) IT(x, t) 



7 (t)|A i (x,*)) = |4CM)+&(x)) 



(13) 



(14) 



It is clear that 7 does not change E and, since H = curl A, it does not change H either. 
Therefore, by gauge-invariance, 



<0| T (*)|0> = <0|0> = 1 
Performing the Wick rotation to Euclidean space, we obtain: 

// E (x, t) = exp I -5 J d 3 y ^(x, y) 11; (x, x 4 ) 



(15) 



(16) 



From now on, we will be interested only in the Euclidean quantity and drop for simplicity 
the E subscript. 

For free photons, rescaling the fields by a factor l/y/J3 with (3 = 1/e 2 , 



IF f + Fnb 



4 r [IV r \iV 



The integral is Gaussian and can be directly computed giving: 

(/*) = exp / ^{Foi{k)F Qj {k))h^) bj(k) 



(17) 



(18) 



Since 
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(F 0i (k)F oj (k)) = M^±M (19) 



if b is such that kibi(k) = 0, we have 



The same calculation for a pure gauge transformation yields: 



^) = ^iU^rMk)\ 2 ) (21) 



2 J (2ny 

We realize that the first term in the exponent of eq.(^) is a normalization which can be 
subtracted by taking instead of (/i) the ratio: 

W = g (22) 
where (7) is defined by means of any gauge transformation g, such that: 

'^)| 2 = /44|b(*0| 2 (23) 



(2vr) 4 ' ov " 7 (2tt 



Then 

(A) = ex P |-^ J^|k||b(*)| a j (24) 

The integral in the exponent, once regularized at small distances, tends to +00 as V — > 00. 
In the infinite volume limit (p) = as it should be, since the perturbative vacuum has 
zero magnetic charge. The same situation appears when computing the overlap of the Fock 
vacuum to the Bogolubov rotated vacuum in a superconductor J23 . 



The same result is obtained if we work directly in the Schrodinger representation in the 
Coulomb gauge. There the vacuum wave functional is: 

tf(fl) = []C k e-^ k l^(k)l 2 (25) 

k 

and 

(n|A*|n> = exp j-i/J J ^|k| |b(A;)| 2 | (26) 
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III. LATTICE FORMULATION 



A lattice version of the operator /j, is obtained by replacing eFoi by the plaquette Po«, or 
better, by its imaginary part: 



eF ( 



Oi 



ImP 0i 



(27) 



and discretizing the field bi. 

Then the disorder parameter becomes: 



VU exp{-[3S} exp {-±f3J2n,ibi(n)Im P 0i (n)} 
JvUexp{-(3S} 

or, if we want to cancel the unwanted normalization, we can divide by: 

VU exp{-/3S} exp{-±PY^i9i(n)I™Poi(n)} 



(28) 



(7) 



obtaining 



JvUexp{-(3S} 
JvU exp {-PS} exp { - 1/3 £ n>i h (n)Im P 0i (n) } 



(29) 



(30) 



JVU exp { -f3S} exp { - \f3 E n ,i 9i {n)ImP 0l {n)} 

We stress once again that any gauge function <?j is acceptable, provided the normalization 
condition (23) is satisfied. 

If we blindly compute (fi) or (p) by numerical simulations, a first technical difficulty 
arises. We are faced with the usual problems encountered in computing quantities like 
a partition function, which are exponentials of extensive quantities, proportional to the 
number of degrees of freedom. The distribution of the values is not Gaussian and the error 
does not decrease by increasing statistics (see, e.g. |24], where the same problem appears in 
a different context). To avoid that, we will compute the quantity: 



P=^log(^) = ^log( / .)-A log(7) 



At (5 — 0, (/i) = (7) = 1, and therefore: 

log(/2) 

We easily get, putting 



n =0 



S 9 = oE9'( n ) /mP ( 



Ot 



"0=0 



(31) 

(32) 

(33) 
(34) 



the following expression for p: 



p= (S + S g )s+s g - (S + S b ) s+Sb (35) 

which can be evaluated by numerical simulations. The subscript on the average indicates 
the action defining the Feynman integral. 

The two quantities on the rhs have the same strong coupling expansion. Thus, the 
use of (p) instead of (//), besides producing a cancellation of the spurious normalization 
of the Feynman path integral, can help in eliminating the lattice artefacts produced by the 
discretization which can spoil the continuum limit. This brings us to a second, more physical 
difficulty, which is the continuum limit. In fact, while for QCD, which is asymptotically free, 
we expect that, at sufficiently high ft, lattice artefacts should cancel, for a model like U(l) 



this point is not so clear in principle. In Ref. [21] a proof is given of monopole condensation 



in the confined phase of U(l), defining a disorder variable for the Villain action. We have 
checked that our (p) operator exactly coincides with the one of Ref. [J21|, when we use the 



Villain action: we expect that for Wilson action the same will hold. We have then computed 
numerically p. 

For a good order parameter, we would expect p to be zero, or (p) = 1 below the critical 
value of the coupling ft c and then to show a large negative peak around ft c , corresponding to 
a drop to zero of (p). At larger values of ft, we have free photons and Eq. fl2"4| ) should hold. 

Figure 1 shows the behaviour of p for a 12 4 lattice, for a monopole in the center of the 
space lattice. In order to be able to identify the signal as a genuine physical result (i.e. not 
due to lattice artefacts), we have performed a number of checks: 

1) we have changed the form of 6j by a gauge transformation to get the Wu-Yang expres- 
sion of the monopole potential. The result does not change qualitatively. 

2) (7) shows practically no signal at ft c within the errors, and does not change appreciably 
by changes of g^. To test that we have computed (7) for two different choices of 
<7i(x), both satisfying the condition Eq. (p3|) : (71) for gi(x) = const, and (72) for 
<72(x) oc x/|x| 2 . In Fig. 2 we display 

- Ai i^l 

Pgauge- ^ ^ 

and p = In for a 6 4 lattice. p gaU ge shows no relevant signal at ft c . 

3) We have measured the correlation function between a monopole antimonopole pair at 
large time distance, in the same position in space. To do that we define 

C(d) = (/x(0,0)/x(0,d)) 



1 



S(b,b)( d ) = -J2 b i( n ) ImP 0i( n ) „ _ -J2 b i( n ) ImP 0i( n 



1 



no=0 



2 



no=d 



and we measure 
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P(b,b)(d) - 2(5 + Sg)s+s g - (S + 5 {6i 5)(d))s+s (il>5) (d) 

Since 

Jp ln J^ = P(b,~b)( d ) c (d)\/3=o = 1 c (d) = J Q P(b,b){d)d/3 

By the cluster property we should have at large d that C(d) — > (p(0, 0)) 2 , or p, h ^.{d) — >■ 
2p. Fig. 3 shows that this expectation is indeed verified. 

We notice that the height of the negative peak of p at /3 C increases with volume [see Fig.l 
and Fig. 2]. The value of f3 c as defined by the position of our peak is f3 c = 1.01(1) for a 6 4 
lattice and (3 C = 1.009(1) for a 12 4 lattice. 

We have taken monopole charge n = 4 (Eq.(|2])) to get a good signal with a relatively low 
statistics (tipically 10 4 configurations per value of 0)\ smaller charges (n = 1,2) give the 
same results but the signals are smaller and more noisy. Our statistical errors are shown in 
the figures, when they are larger than the symbols used. 



IV. CONCLUSIONS 

The strong negative peak in p at f3 c is a clear signal that monopoles are condensed in 
the confined phase, while they are not in the high-/3 regime of the U(l) model. The same 
signal can be extracted from (p p) correlators. We expect a similar signal in QCD at the 
deconfining phase transition for the disorder variable creating U(l) monopoles defined by 
abelian projection. Work is in progress and a comprehensive report will appear soon. 

We underline that the technique used is exactly the same as for the U(l) and that it can 
be extended straightforwardly to different systems in which the condensation of solitons can 
play a physical role, like the 3-d XY model, or the 4-d Georgi-Glashow model. Work is in 
progress also on these systems. 
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FIGURES 
FIG. 1. p versus (3 on a 12 4 lattice 

FIG. 2. p and p gaU ge versus (3 on a 6 4 lattice. 

FIG. 3. /^(cf) versus (3 at distances d = 4, 7, 9, compared to 2p. 
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